Ultracold Rydberg Atoms in a loffe-Pritchard Trap 
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We discuss the properties of ultracold Rydberg atoms in a loffe-Pritchard magnetic field configu- 
ration. The derived two-body Hamiltonian unveils how the large size of Rydberg atoms affects their 
coupling to the inhomogeneous magnetic field. The properties of the compound electronic and cen- 
ter of mass quantum states are thoroughly analyzed. We find very tight confinement of the center 
of mass motion in two dimensions to be achievable while barely changing the electronic structure 
compared to the field free case. This paves the way for generating a one- dimensional ultracold 
quantum Rydberg gas. 

PACS numbers: 32.60.-|-i,33.55.Be,32.10.Dk,33.80.Ps 



I. INTRODUCTION 

Powerful experimental cooling techniques have been 
developed in the past decades that allow us to probe the 
micro and nanokelvin regime while controlling the inter- 
nal and external degrees of freedom of atomic systems. 
As a result dilute ultracold gases that qualify perfectly for 
the study of quantum phenomena on a macroscopic scale 
P, 0, Q can nowadays be prepared almost routinely. Al- 
though being dilute, interactions play an important role 
and rich collective phenomena, reminiscent of e.g. those 
in traditional condensed matter physics, appear. 

The attractiveness of Rydberg atoms arises from their 
extraordinary properties [j]. The large displacement of 
the valence electron and the atomic core is responsible 
for the exaggerated response to external fields and, there- 
with, for their enormous polarizability. Rydberg atoms 
possess large dipole moments and, despite being electron- 
ically highly excited, they can possess lifetimes of the 
order of milliseconds or even more. 

Due to their susceptibility with respect to external 
fields and/or their long range interaction, ensembles of 
Rydberg atoms represent intriguing many-body systems 
with rich excitations and decay channels. Starting from 
laser cooled ground state atoms, a laser typically excites a 
subensemble of the atoms to the desired Rydberg states. 
Since the ultraslow motion of the atoms can be ignored 
on short timescales, Rydberg-Rydberg interactions domi- 
nate the system and we encounter a so-called frozen Ryd- 
berg gas [5|. The strength of the interaction can be var- 
ied by tuning external fields and/or by selecting specific 
atomic states. An exciting objective arc the many-body 
effects to be unraveled in ultracold Rydberg gases (see 
Refs. 0] and references therein). At a certain stage 
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of the evolution ionization might take over leading to a 
cold Rydberg plasma. 

Beyond the above there is a number of topical and 
promising research activities involving cold Rydberg 
states. One example are long range molecular Rydberg 
states |8| with unusual properties if exposed to magnetic 
fields Q. Another one is due to the strong dipole-dipole 
interaction of Rydberg atoms which strongly inhibits ex- 
citation of their neighbors [HI ■ The resulting local 
excitation blockade is state dependent and can turn Ryd- 
berg atoms into possible candidates for quantum infor- 
mation processing schemes [l^, [13] ■ 

A precondition for enabling the processing of Rydberg 
atoms is the availability of tools to control their quantum 
behavior and properties. An essential ingredient in this 
respect is the trapping of electronically highly excited 
atoms. The present work provides a major contribu- 
tion on this score. Let us briefiy address previous works 
on Rydberg atoms exposed to inhomogeneous static field 
configurations. 

First evidence for trapped Rydber g g ases has been 
experimentally found by Choi et al. |l4l . [Tsj . The au- 
thors use strong bias fields to trap "guiding center" drift 
atoms for up to 200 ms. Quantum mechanical studies 
of highly excited atoms in magnetic quadrupole fields 
demonstrated the existence of e.g. intriguing spin po- 
larization patterns and magnetic field-induced electric 
dipole moments [H, [l3| • These investigations were based 
on the assumption of an infinitely heavy nucleus. A 
description of the coupled center of mass (c.ni.) and 
electronic dynamics has been presented in Refs.p^. [l9|: 
Trapping has been achieved for quantum states with suf- 
ficiently large total, i.e. electronic and cm., angular mo- 
mentum. Pictorially speaking this addresses atoms that 
circle around the point of zero field at a sufficiently large 
distance. Recently it has been demonstrated that trap- 
ping in a loffe-Pritchard configuration is possible without 
imposing; the condition of large cm. angular momenta 
[20| . The present investigation works out this setup in 
detail and provides comprehensive results for Rydberg 
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atoms exposed to the loffe-Pritchard field configuration. 

In detail we proceed as follows. Sect. |TT] contains a 
derivation of our working Hamiltonian for a highly ex- 
cited atom in the inhomogcncous field including the cou- 
pling of the electronic and cm. motion of the atom. In 
Sect, mil we introduce an adiabatic approximation in or- 
der to solve the corresponding stationary Schrodinger 
equation. In Sect. IIVI we analyze the obtained spectra 
and point out the capacity of the loffe bias field to regu- 
late the distance between the surfaces, and with that the 
quality of the adiabatic approach. Intersections through 
the surfaces show their deformation when the field gradi- 
ent is increased. Subsequently we characterize the elec- 
tronic wave functions by discussing relevant expectation 
values. Sect. |V] is dedicated to the cm. dynamics in 
the uppermost adiabatic energy surface. We arrive at 
a confined quantized cm. motion without the need to 
impose any restriction on its properties. Examining the 
fully quantized states we observe that the extension of 
the electronic cloud can exceed the extension of the cm. 
wave function. 



II. HAMILTONIAN 
A. Two-body Approach 

The large distance of the highly excited valence elec- 
tron (particle 1) from the remaining closed-shell ionic 
core of an alkali Rydberg atom (particle 2) renders it 
possible to model the mutual interaction by an effective 
potential which is assumed to depend only on the dis- 
tance of the two particles. For alkali atoms, in partic- 
ular, whose core possess zero total angular momentum 
and zero total spin, the only essential difference to the 
Coulombic case is due to the finite size of the core. In any 
case, the effective potential V{r) only noticeably differs 
from the pure Coulomb potential at small distances r. 
States of high electronic angular momenta I, on which 
we focus in the present investigation, almost exclusively 
probe the Coulombic tail of this potential. 

The coupling of the charged particles to the external 
magnetic field is introduced via the minimal coupling, 
p ^ p — qA, where q is the charge of the particle and A 
is a vector potential belonging to the magnetic field B. 
Including the coupling of the magnetic moments to the 
external field (fii and /X2 originate from the electronic 
and nuclear spin, respectively), our initial Hamiltonian 
reads (we use atomic units except when stated otherwise) 

H^n^t = ^ (pi - 9l A(ri ))' + ^ (p2 - <Z2 ^(rs))^ 

- r2\) - /XI • B(n) - fi2 ■ B{r2) . (1) 

We do not take into account spin-orbit-coupling and rel- 
ativistic mass changes. The difference in energy shift 
for adjacent, large angular momentum states {I, I ± 1) 
due to these relativistic corrections is /SWfs = Q!^/2ri'^ 



j2ll |. where a is the fine structure constant, and therefore 
negligible for Rydberg states. At n = 30 one receives 
l^WpS = 1.1 X 10"^^ atomic units. To give an idea 
of the scope of this approximation we anticipate a re- 
sult from Sec. IIVI The energy gap between two adjacent 
high-1 electronic states is approximately Edist = B /2 a.u. 
Demanding ISWps /Edist *C 1 results is constraining the 
loffc field strength B to be much larger than 5 niG. 

Before we focus on the loffe-Pritchard configuration let 
us first examine a general field B composed of a constant 
term Be, a linear term Bi and higher order terms, B = 
^ Bi . The vector potential shall satisfy the Coulomb 
gauge. The squared terms can then be simplified tak- 
ing advantage of the vanishing commutator [A(ri),pi] 
to obtain (pi — qA{ri)Y = p\~ '^qA{ri) pi + q'^ A{ri)'^ . 
In the so-called symmetric gauge the vector potential of a 
constant magnetic field is given by Ac(ri) ~ 1/2 Be x ri. 
The analogon for a linear field is A; (ri ) = l/3JB((ri)xri. 
It can be proven that the vector potential of an arbitrary 
magnetic field can be expanded in a corresponding form 
[2^ permitting a representation of the vector potential as 

a cross product A{ri) — X^i^K^i) = -^(^1) ^ '''i^ where 
B{ri) = "^giBilri) and i e {c, Z,...} denotes the or- 
der of the corresponding terms of A and B with respect 
to spacial coordinates, gi are the coefficients i, ^ etc. 
The particular form of this potential and the vanishing 
divergence of magnetic fields admit the simplification 

A{ri)-pi = {rixpi)-B{ri)=Li-B{ri) , (2) 

where we cxemplarily defined the angular momentum of 
particle 1, Li = ri x pi. 

Since the interaction potential depends only on the 
distance of the two particles, it is natural to introduce 
relative and cm. coordinates, ri = /? + {M2/M)r and 
r2 = R-{Mi/M)r with the total mass M = Mi+Kh- If 
no external field was present, the new coordinates would 
decouple the internal degrees of freedom from the exter- 
nal cm. ones. Yet even a homogeneous magnetic field 
couples the relative and the cm. motion [23, [23|. For 
neutral systems in static homogeneous magnetic fields, 
however, a so-called 'pseudoseparation' can be performed 
providing us with an effective Hamiltonian for the rela- 
tive motion, that depends on the cm. motion only para- 
metrically via the eigenvalues of the pseudomomentum 
0, [2^, m, [13 which is associated with the cm. motion. 
Such a procedure is not available in the present case of a 
more general inhomogeneous field. In the new coordinate 
system the Hamiltonian ([1]) becomes 

H = Ho + LiB{R + ^r) - L^BiR - ^r) 
-p,^B{R+^r)-p,2BiR-^r)+OiA') , (3) 

where the angular momenta of the particles read 

Li = {Mi/M)Lii + {M2/M)Lr + Rxp+ {m/M)r x P 
L2 = {M2/M)Lr + {Mi/M)Lr - Rxp {m/M)r x P 
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(see also Ref. [3), and the terms that do not depend 

on the field are summarized to Hq = + ^ + V{r). 
Here, L^. = r x p, Lji = R x P, and the reduced mass 
m = M1M2/M have been introduced. 

To simplify the Hamiltonian we apply a unitary trans- 
formation that eliminates cm. momentum dependent 
coupling terms generated by the homogeneous field com- 
ponent 



U = exp < - Be X r ■ R 



(4) 



Ha transforms as follows 



U^HoU = H0 + -B, 



-Rx p 



M 



-r X P 



-OiBf. 



The transformation of the remaining terms generates ex- 
clusively additional terms, that are quadratic with re- 
spect to the magnetic field. Exploiting now the fact that 
the mass of the ionic core is much larger than the mass of 
the valence electron, we only keep magnetic field depen- 
dent terms of the order of the inverse light mass 1 /Mi 
(which becomes 1 in atomic units). We arrive at the 
Hamiltonian 

U'^HU = p^/2 + C/V(r)[/ + P'^/2M + 1/2 L^- Be 
+ Ai{R + r)-p+iLr+Rxp)-Bn{R + r) 

- fii - B{R + r) fj,2 B{R) . (5) 

The diamagnetic terms which are proportional to 
(and herewith proportional to B^, see Eq. ([2])) have 
been neglected. Due to the unitary transformation U, 
i2-dependent terms that are quadratic in the loffe field 
strength B do not occur and only an electronic term 
B^{x^ + y^)/8 remains whose typical energy contribu- 
tion amounts to B^n^/8 « lO^B^ for n = 30. Besides 
we obtain a term quadratic in the field gradient G. The 
term quadratic in the loffe field is negligible in compar- 
ison with the dominant shift due to the linear Zeeman 
term as long as B is significantly smaller than 10^ Gauss 
which is guaranteed in our case. Moreover, the c. m. 
coordinate dependance of this diamagnetic term is much 
weaker than the c. m. coordinate dependance of the terms 
linear in the field gradient. The term quadratic in the 
field gradient can be neglected in comparison with the 
corresponding linear term. Up to now we did not use 
the explicit form of the loffe-Pritchard field configura- 
tion. (In anticipation of the special field configuration 
we leave the term containing Ai in its original form.) 



B. loffe-Pritchard Field Configuration 

Two widely spread magnetic field configurations that 
exhibit a local field minimum and serve as key ingre- 
dients for the trapping of weak-field seeking atoms are 
the 3D quadrupolc and the loffc-Pritchard configuration. 



The loffe-Pritchard configuration resolves the problem of 
particle loss due to spin flip by means of an additional 
constant magnetic field. A macroscopic realization uses 
four parallel current carrying loffe bars which generate 
the quadrupolc field. Encompassing Helmholtz coils cre- 
ate the additional constant field. There arc many alter- 
native layouts, the field of a clovcr-lcaf trap for example 
features the same expansion around the origin [2^. On 
a microscopic scale the loffe-Pritchard trap has been im- 
plemented on atom chips by a Z-shaped wire poj . 
The vector potential and the magnetic field read 



(6) 




(7) 



where A„ 



-I- xey) and Bq = 
Br is the constant 



Q{2xzex + 2yzey + {x^ + y^ — 2z'^)ez). 
field created by the Helmholtz coils with B being the 
loffe field strength. Bi originates from the loffe bars 
and depends on the field gradient G. Bq designates the 
quadratic term generated by the Helmholtz coils whose 
magnitude, compared to the first Helmholtz term, can 
be varied by changing the geometry of the trap, Q = 
B ■ |(i?2 _ 4Li2)/(i?2 ^ £)2)2^ .^j^gj.g ^ ig ^j^g x&dins of 

the Helmholtz coils, and 2D is their distance from each 
other. 

If we now insert the special loffc-Pritchard field con- 
figuration using Eqs. (j6l7p into the transformed Hamil- 
tonian ([5|) we obtain 

Hip = Ha + P^ /2M + BLJ2 + G{x + X){y + Y)p^ 
+ Q/^[{{x + X)py-{y + Y)p.,) 

■ {{x + Xf + {y + Y - 2{z + Z)){y + Y + 2{z + Z)))] 
- tiiB{R + r)- ti2BiR) , (8) 

where Ha = — 1/r is the operator for a field free 
atom. The well known Zeeman term BLz/2 comes from 
the uniform loffe field generated by the Helmholtz coils. 
The following term, involving the field gradient G, arises 
from the linear field generated by the loffe bars and 
couples the relative and cm. dynamics. The part in 
the squared brackets originates from the quadratic term, 
again created by the coils. It is the only one that de- 
pends on the Z coordinate, we will see below that its 
contribution is negligible under certain conditions. The 
last term couples the spin of particle two to the magnetic 
field. Since the electronic spins of closed shells combine 
to zero, the spin of particle two is the nuclear spin only. 
Even though scales with I/M2, we will still keep 

the term. Being the only one containing the nuclear spin 
it is essential for a proper symmetry analysis. 
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Operator 




operation 




Px 


X parity 


X — > — X, }C 


-> -X 




electronic spin x op. 


Sy > Sy , 




V 


nuclear spin x op. 


Ey — > — Ej,, 


V i V 


Ixy 


coordinate exchange 


X y, X ^ 




Sxy 


el. spin component exc. 


Sx ^ Sy , 


Sy > Sx 




nuclear spin comp. exc. 


Ex ^ ^y: 




T 


conventional time reversal A ^ A* 





TABLE I: Symmetry operation nomenclature. ± j, uj, emu 



are exemplified by j 



Pj , Sj , and Ej 
X, but hold of course also for j = y, z. 



C. Symmetries, Scaling and the Approximation of 
a Single n-Manifold 

Our Hamiltonian is invariant under a number of sym- 
metry transformations Us that are composed of the el- 
ementary operations listed in Tab. [H The parity oper- 
ations Pj, j £ {x,y,z}, are defined by their action on 
the spatial laboratory coordinates of the particles which 
translates one-to-one to cm. and relative coordinates. 
In order to exchange the x and y components of the elec- 
tronic spin we introduce the operator 



Sxy — 



where SxyS*y = 1. T represents the conventional time 
reversal operator for spinless particles which, in the spa- 
tial representation, corresponds to complex conjugation. 
Our unitary symmetries are 



PxPySz^z 

PyPz-^xy^xy^xy 

PPT q* Y* 

^xJ^ z^xyi-'xy'^xy ■ 



(9a) 
(9b) 
(9c) 



The Hamiltonian is also left invariant under the antiuni- 
tary symmetry transformation 



TR, 



(10) 



By consecutively applying the latter operator and the 
unitary operators ([5a|) . (j9bp and ([5c|) it is possible to 
create further antiunitary symmetries: 



TPxSz^z 

PPz ^xySxy^xy 

TP P P J V* 

^ ^x-Tyr z-LxyJxy^xy 



(11a) 
(lib) 
(11c) 



and T?^y 



Paying regard to the fact that S^y = —Sz 
— and that T neither commutes with 5*1, nor with Sxy 
and 'Sxy, one finds that the operators (|9aj|llcp form a 
symmetry group. 

If no loffe field is present {B = 0), eight additional 
symmetries can be found leaving the Hamiltonian invari- 
ant. For an effective one particle approach (and the cor- 
responding one particle symmetries) this was discussed 
in Ref. 13011. 



As indicated before, the quadratic magnetic field term 
is small and can be tuned by changing the trap geometry. 
It can provide a longitudinal confinement which may be 
treated by perturbative methods. In the case of negligi- 
ble quadratic field Bq, which we assume in the following, 
the term in the squared brackets of the Hamiltonian Q 
drops out and the Z coordinate is cyclic. The correspond- 
ing conjugated momentum Pz is consequently conserved 
and the longitudinal motion is integrated by simply em- 
ploying plane waves \kz) = exp{iZkz}- The constraints 
for this approximation to be valid can be obtained by 
comparing the above-mentioned term in squared brack- 
ets with the Zeeman term, BLz/2. Estimating {x) ~ , 
(xpy) w (ypx) ~ n, and using \Q\ £ B/{D'^ + R^) we find 



L*^ + i?^ > and 



^n(D2 + i?2) » X, r , 



(12) 
(13) 



where D and R characterize the trap geometry. 
Eqs. (jl2ll3p are easily fulfilled. We are therefore left 
with the Hamiltonian 



H = Ha + (P^ + P^)/2 + He , 
where the electronic Hamiltonian reads 



(14) 



H, = BLz/2 + G{x + X){y + Y)pz-l^iB{R + r) . (15) 

For all laboratory fields one finds the magnetic field 
strength B and the magnetic field gradient G to be a 
lot smaller than 1. Our Hamiltonian ([8]) is thus dom- 
inated by Ha- The energies of the field free spectrum 
E\ = —1/2^? are n^-fold degenerate. We can assume the 
loffe-Pritchard field not to couple adjacent n-manifolds 
as long as \E\ - E'^'^\/Ezee > 1- The resulting con- 
straints B < 71"**, G < 71"'^ and GR < n^'^ yield 
B < 2900 G, G < 6 • 10^ T/m for ?i = 30 and 
i? <C 2.9 mm if we additionally assume the field gradient 
G to be as large as 100 T/m. In our parameter regime 
each 7T.-manifold can therefore be considered separately. 
We thus project the full Hamiltonian on the hydrogenic 
eigenfunctions \a) = \n,l,mi,ms), HaIo) = E'^\a), with 
fixed principal quantum number n, that cover an entire n- 
manifold. I denotes the orbital angular momentum quan- 
tum number, m; the one of its z component Lz and rris 
stands for the quantum number of the electronic spin. 

Working in a single n-manifold we can reformulate the 
term in the Hamiltonian (|14p involving the field gradi- 
ent G into a more compact form. We first consider the 
commutator [yz,HA] = [yz,p'^]/2 = i{ypz + zpy). This 
yields 



{a\ypz\a 



{a\zpy\a') = -i{a\[yz,HA]\a.') = , (16) 



since |a) and \a') are eigenkets to the same eigenvalue 
En- Establishing the relation to the orbital angular mo- 
mentum operator via ypz = Lx + zpy results in 



1 



{{a\ypz\a')) ^ -{{a\Lx\a')) . 



(17) 
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The same procedure can be applied to xpz leading to 



{{a\xpAa')) = -\{{a\Ly\a')) 



(18) 



Furthermore {a\XYpz\a') = since p, ~ and 
eventually wc can write 

G{x + X){y + Y)p, = G{xyp, + XL^/2~YLy/2) , (19) 

where we omitted the bracketing alphas, but keep in mind 
that the above identity holds in a single n-manifold only. 

In order to remove the separate dependencies on the 
field parameters B, G, and on the mass M from the 
coupling terms, we introduce scaled cm. coordinates, 
R 7^ 3 R, with 7 = GM, and simultaneously we intro- 
duce the energy unit 6 = 73 /M. Introducing the effective 
magnetic field 



G{X,Y) 




(20) 



and omitting the constant energy offset £"2, the Hamil- 
tonian can be given the advantageous form 



n 



p2 

V 



-fi-G{X,Y)+T'ixyp,+xS,~ySy). (21) 



The first term is the cm. kinetic energy, is the 2n^- 
dimensional matrix representation of the total magnetic 
moment of the electron, ^{Lr + 2S), and the second term 
in ([2T|) describes its coupling to the effective magnetic 
field G. The latter results from the original field Bc + Bi 
in Eq. ([7]) taking into account the corresponding coordi- 
nate and energy scaling factors. Si are the components 
of the electronic spin, S = —fii. The nuclear spin term 
— /X2 • B{R) has been omitted since it is several orders of 
magnitude smaller than the electronic one. 



III. ADIABATIC APPROACH 

The large difference of the particles' masses and veloc- 
ities in our two body system makes it plausible to adi- 
abatically separate the electronic and the cm. motion. 
The corresponding time scales differ substantially even 
for large principal quantum numbers n. However, due 
to the enormous level density in case of Rydberg atoms 
it is a priori unclear whether isolated energy surfaces 
might exist or whether, as one might naturally assume, 
non-adiabatic couplings are ubiquitous and therefore an 
adiabatic approach might invalidate itself. The proce- 
dure is reminiscent of the Born-Oppenheimer ansatz in 
molecular systems and is based on the idea that the slow 
change of the heavy particle's position allows the elec- 
tron to adapt instantaneously to the inhomogeneous field. 
The electronic energy of the system can thus be consid- 
ered as a function of the position of the heavy particle. 



The adiabatic approximation is introduced by sub- 
tracting the transversal cm. kinetic energy, T — 
(P2 + P^)/2, from the total Hamiltonian The re- 

maining electronic Hamiltonian for fixed center of mass 
reads 



Tie ^ 11- G{X, Y) + 73 {xyp^ + xSa; ~ ySy 



(22) 



The electronic wave function ip,^ depends parametrically 
on R and the total atomic wavefunction can be written 
as 

\^{r,R)) = \^,{r;R)}^\MR)) , (23) 

where \^i,{R)) is the center of mass wave function. 
The internal problem posed by the stationary, electronic 
Schrodinger equation 

He \^.{r; R)) = E,{X, Y) |^„(r; R)) (24) 

is solved for the adiabatic electronic potential energy sur- 
faces Ek,{X, Y), that serve as a potential for the c. m. dy- 
namics. Within this approximation, the equation of mo- 
tion for the center of mass wave function reads 



{T + Ejx,Y)) \MR)) = ^u \MR)) 



(25) 



The spatially dependent transformation U {X, F), that 
diagonalizes the matrix representation Tie of the elec- 
tronic Hamiltonian, is composed of the vector represen- 
tations of the electronic eigenfunctions, = {Uko) = 
{{a\ifi^{r] R))). Since Z// depends on the cm. coordinates, 
the transformed kinetic energy involves non-adiabatic 
couplings AT 

U^nU = U^njA +U^TU =^ E^{X,Y) + T + AT (26) 

that have been neglected in the adiabatic approximation 
of Eq. 



(27) 



AT ^ -1/2- {U^idjiU) +11^8^11) 

+ 2U\dxU)dx + 2U\dYU)dY) 



They can be calculated explicitly as soon as the electronic 
adiabatic eigenfunctions have been computed. Non- 
adiabatic contributions can be neglected if the conditions 



^ {ip^,\{dx'H)\ip,, 
Em — Ek, 

Ek.1 — Ek, 



< 1 



< 1 



(y>.,|(9yH)|^.) 

(^.,|(9f.7^)|^.) 
Em — Ek 



l«l, (28) 
I < 1 (29) 



are fulfilled [19|. The energy denominator in ([28)l and 
(|29p indicates that one can expect non-adiabatic cou- 
plings to become relevant between the adiabatic energy 
surfaces when they come very close in energy, i.e. in the 
vicinity of avoided crossings. 

Recalling the results of the symmetry analysis, it 
can be demonstrated that the energy surfaces E^, ex- 
hibit three mirror symmetries. Within the adiabatic ap- 
proximation, X and Y are parameters in the electronic 



6 



Schrodinger equation. Symmetry operations applied to 
the electronic Hamiltonian thereby merely act onto the 
electronic subspace. If we apply the corresponding re- 
stricted symmetry operation Up = PxPySz^z (Pa|) . that 
was already shown to leave the full loffe-Pritchard Hamil- 
tonian IH]) invariant, to the electronic Hamiltonian He 
(fTS]) . we find 

Ul,Heir;X,Y)Up ^ Heir;-X,-Y) . (30) 

Since unitarily equivalent observables, A and WAU, 
possess the same eigenvalue spectrum, we find the en- 
ergy surfaces to be inversion symmetric with respect to 
the origin in the X-Y plane. The symmetry operator 
Uy ~ TPy, and the operator that is composed of Uy and 
Up, namely Ux = TP^Sz^z (see ^ and pla|l l. mirror 
the energy surfaces at the axes, 

UlHe{r;X,Y)UY = He{r;X,-Y) , (31) 
Uj,Heir;X,Y)Ux ^ Heir;-X,Y) . (32) 

The electronic problem (|24p . with the core fixed at an 
arbitrary position, is three-dimensional. No symmetry 
arguments can be exploited to reduce the dimensionality 
of the problem. In order to solve it, we employ the vari- 
ational method, which maps the stationary Schrodinger 
equation onto an ordinary algebraic eigenvalue problem. 
Since the matrix representation of the electronic Hamil- 
tonians is sparsely occupied, an Arnoldi decomposition 
is used. Both, this decomposition and the surfaces' mir- 
ror symmetries, help to reduce the computational cost of 
solving the electronic Schrodinger equation. 



IV. ELECTRONIC POTENTIAL ENERGY 
SURFACES 

In this section the properties of the electronic adiabatic 
energy surfaces are analyzed for different regimes of loffe 
field strengths and field gradients. These two parameters 
can be used to shape the potential in which the center 
of mass dynamics takes place. To understand how this 
takes place, we inspect the electronic Hamiltonian to un- 
ravel the influence of the individual terms for different 
parameter regimes. 

The characteristic length scale of the center of mass 
dynamics is of the order of one in scaled atomic units. 
It is therefore adequate to compare the magnitudes of 
the different parts of the electronic Hamiltonian ([22]) in 
order to estimate their impact on the center of mass mo- 
tion, putting X and Y equal to one. The first part, 
fi-G{X, Y), consists of the coupling terms X{^Lx + Sx) — 
Y{^Ly + Sy), that are then of the order of {Li) « n for 
high angular momentum states, and of the Zeeman term 
Ci^Lz + Sz), which can be as large as (n. The second 
part, "f^^^{xypz + xS^ — ySy), is quadratic in the rela- 
tive coordinates which makes it particularly important 
for high principal quantum numbers n. If we consider 
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0.064 0.138 0.296 0.638 1.375 2.963 6.383 

0.177 0.382 0.823 1.773 3.820 8.229 17.729 

0.454 0.978 2.107 4.539 9.778 21.067 45.387 
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TABLE II: Explicit values for 7^/^71^ = {GMy/^n'^ and C = 
BM^/^G"^/^ for *''Rb in atomic units. The first block lists 
for different values of the field gradient G and for dif- 
ferent principal quantum numbers n. The second block lists 
for different field gradients and for different field strengths B. 



the expectation values of the relative coordinates to be 
of the order of n^, and {ypz) ~ {Lx) ~ n, the overall 
magnitude can be estimated to In a nutshell, we 

have for the mentioned three terms the following relative 
orders of magnitude. 



1 



and 7 3 n 



(33) 



Due to the special form of the electronic Hamiltonian, 
changing the magnetic field parameters B and G while 
keeping their ratio C/l^^^ = B/G (and n) constant re- 
sults in a mere scaling of the cm. coordinates. We pro- 
vide typical examples for values of the quantities ()33|) in 
table ini 



A. Regulating Capacity of the loffe Field 

To understand the impact of the loffe field strength B 
on the adiabatic energy surfaces, we isolate its effect by 
suppressing other influences. This can be done by choos- 
ing a relatively low field gradient G and/or a small prin- 
cipal quantum number n (sec Tab. llT|) . The factor 7371^ 



becomes small, and the last term in Eq. (|22p will hardly 
provide any contribution. Within this regime, that we 
focus on in this subsection, approximate analytical ex- 
pressions for the electronic adiabatic energy surfaces can 
be derived. We diagonalize the approximate electronic 
Hamiltonian 



He ^\g{L + 2S) . 



(34) 



by applying the spatially dependent unitary transforma- 
tion 



Ud{X,Y) 



g#(L,+S,)g»/3(L^+S„) 



(35) 
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with (f> = arctan y , cos (3 = 
sin/? = -\/X2TT2|G(X,y) 

1 



\M2B\G{X, f)|-i 

This yields 



and 



U],H,Ud = -{L, + 2S,)\G{X,Y) 



(36) 



for the transformed approximate electronic Hamiltonian. 
The spatially dependent transformation Ud locally ro- 
tates the magnetic moment of the electron, which in- 
cludes its spin and its angular momentum, such that it is 
parallel to the local direction of the magnetic field. The 
operators and Sz are not identical to the ones before 
having applied the transformation (j35p . they are rather 
related to the local quantization axis defined by the local 
magnetic field direction [l8l |. 

The adiabatic potential surfaces evaluate to 



-{mi + 2m,)^X^ + + C 



(37) 



The possible combinations of m; and yield 2n + 1 
energy surfaces. The surfaces highest and lowest in en- 
ergy correspond to circular states, {\mi\ = Imax = n — 1, 
mi + 2ms ^ i'^); a-iid they are the only non-degenerate 
ones. For the other surfaces (jm; -I- 2ms\ < n), the 
multiplicity of {mi + 2ms), and with that the degree 
of degeneracy of the corresponding surfaces, is given by 
2n — \mi + 2ms + 1| — + 2ms — 1|- Starting from the 
highest energy surface, the levels of degeneracy thus are 
1,2,4,6,....^ 

The approximate surfaces (|37p arc rotationally 
symmetric around the 0-axis. An expansion around this 
axis {p = VA^ + y^ ^ ^) yields a harmonic potential, 



1 



1 



En{p)~{C+T77P )-n{r'M + 2ms) , 



2C 



while wc find a linear behavior, 



EM 



{ini + 2ms) , 



(38) 



(39) 



when the center of mass is far from the z-axis (p 3> C)- 

For reasons of illustration we demonstrate the behav- 
ior of the adiabatic surfaces with increasing loffe field 
by means of a somewhat artificial example where other, 
previously neglected interactions might be more impor- 
tant. Fig. [1] shows sections through all the surfaces for 
ri = 3. This principal quantum number has been cho- 
sen in order to keep the sections simple while displaying 
the entire n-manifold. We employ ^''Rb in this exposi- 
tory example although the electronic ground state of its 
outermost electron is 5s. The sections have been calcu- 
lated for the field gradient G = 1 T/m and for different 
field strengths B using the total electronic Hamiltonian 
([22)1 . These parameters yield 

^l/3„2 

= 0.003, and values 



for C ranging from 0.01 to 1. The surfaces in the different 
graphs of Fig. [1] indeed validate the approximate expres- 
sion ([57)1 : We find 2ri -I- 1 degenerate surfaces and the 
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FIG. 1: Sections along the X-axis through the electronic adi- 
abatic energy surfaces of an entire n — 3 manifold. The field 
gradient is fixed at G = 1 Tesla/m in order to suppress the 
influence of the last term in He (|22p . From left to right, 
= BM-y~^^^ increases due to an increasing lofi^e fleld. 
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FIG. 2: Sections along the X-axis through the uppermost 
21 surfaces of the n — 30 manifold of ^'^Rb for increasing 
ratios B/{Gn^). The field gradient is fixed at G = 10 T/m 
while the loffe field is increased from top left to bottom right. 
{B = 24 mG, B = 48 mG, B = 0.24 G, B = 0.48 G). For 
small ratios B /{G-n?) the influence of the second term in (|22p 
is not completely suppressed as can be seen from the lifted 
degeneracies in the upper subflgures. 



harmonic behavior for |X| <C C gives way to a linear in- 
crease for \X\ ^ C. The energetic distances and lengths 
in the different graphs are comparable, since the scaling 
factor for the center of mass coordinates 7 = C2M has 
not been changed. We can conclude that increasing the 
loffe field strength B separates the surfaces from each 
other. 
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11.101 
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13.40 
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0.638 
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10 


0.638 


67.006 


0.002 
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1340 


10 


0.638 


670.07 


0.001 


0.01 


6.220 


1 


0.296 


3.107 


0.104 


0.1 


62.20 


1 


0.296 


31.101 


0.001 


1 


622.0 


1 


0.296 


311.022 


0.000 



TABLE III: Minimal distance AE of the two uppermost sur- 
faces of the 71 = 30 manifold. A denotes the discrepancy 
between AE and the approximate predicted value for the dis- 
tance, |, according to Eq. (|40)) . 



The data presented in Fig. [2] have been computed for 
the n ~ 30 manifold. In order to keep the last term in 
small, the field gradient has been set to G = 0.1 T/m 
{—^ ^^/^Ti? = 0.14). The uppermost 21 energy sur- 
faces are shown for different values of the magnetic field 
strength B. Similar to the n = 3 case, one can see the 
harmonic behavior around the origin. The surfaces' min- 
imal distance becomes larger for increasing C,. Since C, 
and are of the same order of magnitude in sub- 

figure (a), the contribution of the last term in ([22|) . that 
lifts the degeneracy of the curves, is visible. 

The energetic distance of the approximate surfaces de- 
scribed by Eq. ([57)1 increases with larger distances from 
the Z^-axis, p, and with larger C,. The minimum energetic 
gap between two adjacent surfaces is at the origin and 
reads 

\E,{0) ~ i?„±i(0)| = |M7-i = ^ . (40) 

The parameter C, (an hence the field strength B) is the 
tool to control the energetic distance between the adi- 
abatic surfaces. Increasing C, one can thus also mini- 
mize the non-adiabatic couplings AT p7|) discussed in 
Sect, mil since they scale with the reciprocal energetic 
distance of the surfaces. 

To check the range of validity of our approximation, 
the minimal energetic distance between the two upper- 
most adiabatic surfaces in the n = 30 manifold has been 
calculated for different parameters, subtracting the full 
2D surfaces from each other, that have been obtained us- 
ing the electronic Hamiltonian (j22[) . One finds the min- 
imal distance to be located at the origin, as expected. 
A in Tab. IIIII denotes the relative deviation between the 
predicted (Eq. (jiU)) ) and the computed value in percent. 
It is small for large loffe field strengths B and low field 
gradients G. Then we have C ^ 7^/'^n^, the last term in 
the electronic Hamiltonian is negligible and our approx- 
imation that leads to is justifiable. 
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X [scaled a.u.] X [scaled a.u.] 

FIG. 3: Sections {Y = 0) through the adiabatic potential 
energy surfaces belonging to the n = 30 manifold of *^Rb 
for decreasing ratios B/{Gv?) = (^/(7^''^n^). The influence 
of the Zeeman term in Tie (|22|) is fixed [C, = 5) while 
increases, (a) B/{G'n) = 10 B = 22.9 mG, G = 4.81 
T/m; (b) B/iGn^) = 5^5 = 91.6 mG, G = 38.5 T/m; (c) 
B/{Gn^) = 1 ^ B = 2.29 G, G = 4807 T/m; (d) draws the 
indicated region in (c) to a larger scale. 

B. High Gradients 

A more complicated picture of the surface properties 
arises when the field gradients become larger. The last 
term in the electronic Hamiltonian, that accounts for fi- 
nite size effects of the atom, 

7^ {xyp;: + xSx - ySy) , (41) 

is no longer small compared to the others in equation 
(P^ . This results in modulations of the adiabatic sur- 
faces we already spotted in the previous section, even 
though the term does not feature any dependency on X 
and Y. These modulations lift the degeneracy that was 
found in the limit of small gradients. Their dependency 
on the cm. coordinates is introduced by the transfor- 
mation U{X, Y) that diagonalizes the electronic problem 
(cf. Secinil. 

In order to isolate the effect of the term (|4T|) on the 
adiabatic surfaces, we vary the scaling factor 7 = GM 
by changing the field gradient G, while keeping ^ = 
BM^/^G^^^^ constant. It is, for example, reasonable to 
demand C, ~ h and to adjust the loffe field strength B to 
meet this condition. Fig. [3] demonstrates the increasing 
influence of the interaction ((1T|) when G is increased. The 
spectra are computed for the n = 30 manifold of ^''Rb, 
C = 5, while G is varied from 4.8 to 4800 T/m. For smafl 
field gradients ((a), B/{Gii?) = 10), the surfaces ap- 
proach the shapes predicted in the limit addressed in the 
previous subsection pV A^ : The adiabatic surfaces with 
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90 100 



FIG. 4: Section through the n = 30 manifold for a field 
strength of 0.01 Gauss and a field gradient of 20 T/m (^''Rb). 
A large number of avoided crossings can be observed. The up- 
permost curve, however, stays isolated from the other curves. 
The insets show the linear behavior of the surfaces far away 
from the 2-axis. 



the same value of the magnetic moment {'mi + 2ms) /2 are 
approximately degenerate. The uppermost energy is the 
only non-degenerate one and to the corresponding eigen- 
state the quantum numbers mi = n — 1 and nis = 1/2 
can be assigned. An increasing field gradient Hfts the 
degeneracy and groups of curves can be observed ((b), 
B/{Gn?) = 5). The energetic distance between these 
groups stays tunable by the bias field strength, as we 
elucidated above (see Eq. (|1D|) '). For even higher field 
gradients, the different parts of the electronic Hamilto- 
nian are of comparable size and finite size effects sub- 
stantially alter the shape of the energy surfaces ((c), 
(d), B/{Gn^) = 1). Avoided level crossings appear and 
non-adiabatic transitions are likely to occur. The up- 
permost energy surface, however, proves to be very ro- 
bust when the field gradient is varied. It is energetically 
well-isolated from the other adiabatic surfaces. Its dis- 
tance to the surface, that is formed by the second high- 
est eigenvalue, only decreases significantly when the ratio 
B/{Gii?) approaches one ((c), (d)). This holds true for 
the entire X-Y-plane. Inspecting the full uppermost sur- 
face one furthermore finds the azimuthal symmetry, that 
is found for large ratios B/{G'n?) (see Sect. IIV A^ . to be 
approximately conserved. 

Another example for the complicated structure of the 
adiabatic electronic energy surfaces is shown in Fig. 31 
The data are calculated for a loffe field strength of 0.01 G 
and a field gradient of 20 T/m. For these parameters, the 
contributions of all terms in the electronic Hamiltonian 
are of the same order of magnitude around X ^ 1. One 
immediately notices the large number of avoided cross- 
ings between the surfaces. The uppermost curve however 
remains isolated from the rest of the curves. Far away 
from the trap center, i.e. for large p = VX^ + Y^, the 



coupling term in X{\Lx + S^) - Y{\Ly + Sy), be- 
comes dominant. A Zeeman like splitting of the surfaces 
emerges, visible in the smaller graphs on the right. 



C. Electronic Wave Functions 

To characterize the electronic wave function ^p^ii"] R), 
that corresponds to the energy eigenvalues constituting 
the uppermost adiabatic surface, we analyze its radial 
extension, angular momentum and spin. The electronic 
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FIG. 5: Expectation value 



of the wave functions that 



correspond to the uppermost electronic energy surface for 
G = 100 T/m (n = 30, *^Rb). B is varied yielding dif- 
ferent values for the ratio C/l^^^n^ = B/Gn^: (a) 0.01 G 
^ B/Gr? = 0.21 a.u., (b) 0.1 G ^ B/Gn^ = 2.1, (c) 1 G 
-> B/Gn^ = 21. The depicted ranges of X and Y correspond 
to 30 characteristic lengths of the cm. motion in scaled units. 



wave function depends parametrically on the cm. posi- 
tion and is, in general, distorted compared to the field free 
case by the external magnetic field. This is reflected in 
the expectation value (r)e(-R) = {'fK.ix;R)\r\ipi^{r\R)) 
which is shown in Fig. [5] for different ratios B/{Gn^). 
The limits of the graphs with respect to X and Y cor- 
respond to thirty characteristic lengths of the cm. mo- 
tion. While keeping G = 100 T/m, B is increased for 
the different plots from left to right. For the smallest ra- 
tio under consideration ((a), B/Gri^ < 1), a pronounced 
maximum of the expectation value (r)e can be observed 
at the trap center. This maximum breaks up into four 
maxima arranged along the diagonals when the ratio is 
increased ((b), B/G-n? > 1), while the amplitude of the 
spatial variation of (r)e decreases. For an even higher 
value of B ((c), BjGr? ^1), only a marginal deviation 
from the hydrogenic field free value for the highest pos- 
sible angular momentum quantum number remains (for 
n = 30 one finds (r)H(n = 30, Z = 29) = 915). In the re- 
gion of local homogeneity, where the magnetic field does 
not vary significantly over the extension of the electronic 
cloud (i.e. far from the z-axis), the expectation value ap- 
proaches the field free value in all subfigures that are 
shown in Fig.O In accordance with the abovementioned 
scaling property of the electronic Hamiltonian Tie, chang- 
ing the field parameters while keeping the ratio BjGr? 
unaltered only modifies the scale of the cm. coordinates, 
whereas the shape of the bright regions and the energy 
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range of the eigenvalues arc not changed. 
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FIG. 6: Expectation values {Lx), (Ly) and (L,) (a,b,c, re- 
spectively) for a ratio B/{G-n) = 2.1 (loffe field B = 0.1 G, 
gradient G = 100 T/m, **^Rb, n=30). In (d) the projection H 
of {Lr) onto the local magnetic field direction G is displayed. 
It is close to the field free maximum value for the angular mo- 
mentum projection, mi^rnax = n — 1. Subplot (e) shows the 
spatial behavior of {L^). The range of X and Y corresponds 
to 30 times the characteristic length of the cm. motion. 

Let us study the angular momentum and its orienta- 
tion. It is to be expected that for dominating lofFe field, 
i.e. for very large ratios B/{Gn'^), the expectation value 
of the angular momentum, (Lr) — {{Lx) , (Ly) , (Lz)), 
is oriented in the loffe-field direction (z-axis). Since 
the loffe field in any case dominates around the origin, 
(Lx) and (Ly) are expected to vanish at {X,Y) = (0,0) 
while (Lz) becomes maximal. This behavior can be 
observed in Fig. [5] where (Li) are displayed (a,b,c) for 
B ~ 0.1 G and G ~ 100 T/m. These parameters yield 
B/{Gn?) = 2.1. The aUgnment of {Lr) and the local 
field direction G{X,Y) is found to be very good in the 
entire X-y-plane (the maximum angle between the two 
is smaller than 3.6°). In subplot (d) we provide the spa- 
tial behavior of the projection of {Lr) onto this local field 
axis, n = {Lr) ■ G{R) /\G{R)\. In the local homogeneity 
limit, n approaches the maximal value for (Lz), namely 
mi^max — n ~ I. In the same manner the expectation 
value (£^), which is displayed in subplot (e), converges 
to the maximal value, lmax{lmax + ^) = n{n—l). Far from 
the z-axis, the uppermost surface hence corresponds to 
the circular state \mi^max, Imax)- The deviation of 11 and 
(L^) from the maximal values close to the z-axis reflect 
the admixture of states with lower quantum numbers m 
and / to the state of the uppermost surface. 

Increasing the applied loffe field by a factor of 10 
{-^ B/{Gv?) = 21), decreases the angle between {Lr) 
and G{X, Y) by a factor of 10^, i.e. a quasi perfect aHgn- 
ment is found. As can be seen in Fig. [71 the projection 
n now only deviates marginally from mi^max- Conse- 
quently, also {L'^) exhibits only minor deviations from 
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FIG. 7: Spatial dependence of the projection 11 of {Lr) onto 
the local field axis for B = 1 G (all other parameters are the 
same as in Fig. [Sjl. For this ratio, B/(Gn^) — 21, the devia- 
tions from the maximal value mi^rnax = n — 1 are marginal. 



(Equally, {L^) 



i.max{lmax + 1), UOt shoWU.) 



its maximum value in the whole X-y-plane. For high 
ratios B/{Gn?), the admixture is therefore marginal and 
one can in a very good approximation assume the elec- 
tronic state in the uppermost surface to be the circular 
state \mi^maxTlmax) for any cm. position. Similar obser- 
vations can be made considering the respective expecta- 
tion values for the spin. For the parameters in Fig. [7] 
the projection of {S) onto G differs less than 10~^ from 
1/2. The expectation values of the examined electronic 
observables converge to the field free values for increasing 
ratios B/{Gn^). 

Our findings indicate that the electronic structure of 
the atom is barely changed in the limit of large ratios 
B/{Gii?). The radiative lifetimes can hence be expected 
to differ only slightly from the field free ones [l^ . 



V. QUANTIZED CENTER OF MASS MOTION 

The energetically uppermost adiabatic electronic en- 
ergy surface is the most appropriate to achieve confine- 
ment. It does not suffer a significant deformation when 
the field gradient is increased and it stays well isolated 
from lower surfaces for a wide range of parameters. Large 
energetic distances to adjacent surfaces suppress nona- 
diabatic couplings (Eqs. (^5]) and (HH])). In order to 
obtain the quantized cm. states we therefore solve the 
Schrodinger equation (|25p for the cm. motion in the up- 
permost surface £'2^2 by discretizing the Hamiltonian on 
a grid. The wave function for the fully quantized state 
is hence composed of the eigenfunction \'Pk{i"i R)) of the 
electronic Hamiltonian in equation (|24p . the wave func- 
tion for the center of mass motion in the X-Y plane, 
|?/'i/(-R)), and the plain wave in Z direction. 



|*(r,H)) = \v.{r;R)) ® |^A.(-R)) <E) \kz) 



(42) 



In Fig. [8] the probability densities of the ground state 
and two excited states of the cm. motion in the upper- 
most surface of the n = 30 manifold of ^^Rb are dis- 
played. These densities reflect the spatial symmetries 
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FIG. 8: Probability densities of the ground state and the first 
and tenth excited states of the cm. motion in the uppermost 
adiabatic potential surface of the n = 30 manifold of *^Rb 
(from left to right). The loffe field strength is set to B = 0.1 
G and the field gradient is G = 10 T/m. 

of the electronic Hamiltonian He ((22)) and consequently 
those of the electronic energy surface. They are com- 
puted for a loffe field strength B = 0.1 G and a field 
gradient of G = 10 T/m, which yields ( = 13.4 and 
C^/^^n? = B/Gn? = 21. According to the discussion 
in Sec. IIV Al the electronic surface then exhibits a har- 
monic behavior around the origin, and the system resem- 
bles the two dimensional isotropic harmonic oscillator in 
the potential Eh{X, ¥) = {(,+ /2C) ■ n/4 (cf. Eq. dM]), 
nil ~ n — \). 

The first two probability densities (from left to right) 
in Fig. [8] explicitely demonstrate the analogy to the har- 
monic oscillator. The nodal structure of the tenth excited 
state is not due to a Cartesian product of ID harmonic 
oscillators but a different combination of the harmonic 
oscillators in the corresponding degenerate subspacc. 

The energies of the cm. wave functions in the approx- 
imative potential Eh{X,Y) read 

= (A^i + A^2 + 1) , iVi,iV2 - 0,1,2... , (43) 

where uP' = n/2C, which are in very good agreement 
with the exact results in the regime where the electronic 
Hamiltonian (|4ip is negligible. Within this approxima- 
tion, the energy level spacing scales with the inverse 
square root of C, Aeh,^ = oj ^ l/\/Ci whereas the ener- 
getic distance of adjacent surfaces scales linearly with 
see Eq. p0|) . 

To describe the properties of the compound quantized 
state, we analyze the extension of the center of mass mo- 
tion, which can be measured by the expectation value 

{p) = {MR)\ ^X^ + y2 i^^(H)) , (44) 

and the mean distance of the core and the electron (r). 
Fig. [9] presents the radial expectation value (p) in Bohr 
radii for the cm. ground state in the uppermost energy 
surface for different parameter sets of the magnetic field. 
For comparison, the expectation value of the cm. state 
in a perfectly harmonic potential, {p)h = -^^o ^ C^^^: is 
also depicted. The char acteris tic length of the cm. mo- 
tion is xq = 1/v^ = {/2Cjn. (Due to the rescaling of 
the cm. coordinates with in Sec. Ill Bl this is of 

the order of 1 for a wide range of parameter sets {B, 
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FIG. 9: Double logarithmic plot of the expectation value 
(p) for the cm. ground state (circles, o) in the uppermost 
adiabatic energy surface (n = 30, *^Rb). The corresponding 
expectation values for the cm. wave function in a perfectly 
harmonic potential are depicted for comparison (+). 



G, n\ in scaled atomic units (cf. Tab. The ex- 

pectation values for the real system, (p), deviate from 
the straight line formed by (p)/i, as the ratio BjG be- 
comes very small. Hence, by choosing large gradients 
and appropriate bias fields, very tightly confining traps 
for highly excited atoms can be obtained (B = 0.1 G aird 
G = 100 T, for instance, give rise to a trap frequency of 
approximately 1.4 MHz). 

The mean distance of the Rydberg electron from the 
core (r) is calculated weighting that very quantity for a 
fixed cm. position (r)e(A', y), with the probability den- 
sity of the cm. wave function: 

(r) ^ {MR)\ {Vn{r;R)\ r |^«(r;H)) \MR)) ■ (45) 

It is depicted in Fig.[Tni along with (p), versus the degree 
of excitation of the cm. motion v. (p) and (r) are of 
comparable size due to the very tight confinement. For 
a loffe field strength of B = 0.1 G and a field gradient 
of G = 100 T/m, for instance, the ratio of (p) and (r) 
for the ground state {i' = 1) is as small as {p)/{r) ~ 0.4. 
The extension of the cm. wave function is thus smaller 
than the extension of the electronic cloud. This strongly 
supports the proposition that our Rydberg atoms cannot 
be considered as point-like particles. 

The expectation value (?') for the electron remains 
nearly constant as the degree of excitation increases, 
and it barely differs from the corresponding field free 
value (dashed line in Fig. pO|) ). As indicated previ- 
ously, we find the electron to be in the circular state with 
nil = n — I, which features the smallest mean square de- 
viation of the nucleus-electron separation (r^) — (r)^ = 
n'^{2n+ l)/4. It is therefore possible, that the cm. and 
the electronic wave function do not even overlap. This is 
indicated in the inset of the upper right plot in Fig. [TU] 
for V = \. 
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FIG. 10: Comparison of the mean extension of the cm. wave 
function, (p), and the mean distance of the core and the elec- 
tron, (r), for the n — 30 manifold of **^Rb. 



VI. CONCLUSION 

We have studied the quantum properties of ultracold 
Rydberg atoms in a lofFe-Pritchard field configuration 
and find trapped cm. quantum states to be readily 
achievable. Our starting point is a two-body approach 
to the Rydberg atom. Relativistic effects and deviations 
of the core potential from the Coulomb-potential as well 
as diamagnetic interactions have not been taken into ac- 
count, which is well justified a posteriori. Applying a 
spatially dependent unitary transformation and addition- 
ally exploiting the major mass difference of the electron 
and the core, we arrived at a two-particle Hamiltonian for 
highly excited atoms in an inhomogeneous field where the 
appearance of the coupling of the relative and cm. dy- 
namics is simplified substantially. Thenceforward we 
have concentrated on the special case of a loffe-Pritchard 
trap. A symmetry analysis of the resulting Hamiltonian 
has been performed revealing seven discrete unitary and 
anti-unitary symmetries. Comparing the energetic con- 
tributions of the different interactions we find it legiti- 
mate to limit our considerations to a single n-manifold 
to solve the corresponding stationary Schrodingcr equa- 
tion. Consequently an adiabatic approach was applied. 
In the ultracold regime the Rydberg electron is much 
faster than the cm. motion of the atom. This justifies 
an adiabatic separation of the internal (relative) and the 
external (cm.) dynamics. The corresponding adiabatic 
electronic potential surfaces have been obtained by diag- 
onalizing the electronic Hamiltonian matrix. In the limit 



of large ratios of loffe field strength and field gradient, 
B/(Gn^), an approximate analytical expression for the 
adiabatic surfaces has been provided. In this limit, the 
surfaces arrange equidistantly and all but the uppermost 
surface are degenerate. The inter-surface distance is then 
proportional to the loffe field strength. The structure 
of the electronic surfaces becomes more complex when 
this ratio decreases. The shape of the uppermost sur- 
face and its energetic separation from others, however, 
prove very robust with respect to changes of the field 
parameters. We hence consider it the most appropri- 
ate to achieve confinement. Exploring the properties of 
the electronic wave functions we find that the expecta- 
tion values approach the field free values when the ratio 
of the field strength and the field gradient, B/{Gn'^), is 
increased. This indicates that, despite the strong local- 
ization of the cm., the electronic structure of the atom 
is barely changed compared to the field free case. Ex- 
amining the compound quantized states we have found 
a regime where the extension of the cm. wave function 
falls below the extension of the electronic cloud, i.e. the 
cm. is stronger localized than the valence electron. In 
this regime Rydberg atoms in inhomogeneous magnetic 
fields can therefore not be considered as point-like parti- 
cles. We conclude that the loffe-Pritchard trap provides 
a strong confinement for Rydberg atoms in two dimen- 
sions that permits their trapping on a microscopic scale. 
For such a one-dimensional guide, a relatively weak lon- 
gitudinal confinement along the z-axis could additionally 
be provided for a non-Helmholtz configuration by the 
quadratic term. As a natural enrichment of the system 
one could study many atoms in that guide. Challenging 
issues are to stabilize such a one-dimensional Rydberg 
gas or to answer the question if it is feasible to use the 
strong Rydberg- Rydberg interaction to create a chain of 
trapped atoms [31| that could then serve as a tool for 
quantum information processing [Tsl . [3^ . [33| making use 
of the state dependent atom-atom interaction. 
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